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Abstract 

We make use of twisted boundary conditions for off-shell, Rome-Southampton renormalisation. 
This allows us to define vertex functions at a fixed physical momentum that need not be a Fourier 
mode. This definition includes choosing a fixed orientation with respect to the lattice axes; only 
then do lattice artefacts, which include 0(4) breaking effects, have a valid expansion in powers 
of the lattice spacing, a. The use of non-exceptional momenta has been found to greatly reduce 
the dependence of the vertex functions on both mass and momentum, p. Excellent statistical 
precision is afforded by plane- wave sources. Together this enables both a theoretically well founded 
and statistically clean continuum limit. Thereafter all dependence can be identified with the 
anomalous running of the operator. We present initial results and develop a practical scheme for 
step-scaling with off-shell renormalization. The size of the steps is continuous rather than discrete, 
allowing arbitrarily small steps and the scheme is easy to implement for general operators . 
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I. INTRODUCTION 



The non-perturbative renormalisation (NPR) of Lattice QCD matrix elements is an im- 
portant topic for the phenomenological relevance of the field. This paper addresses two key 
limitations of the popular Rome- Southampton RI/MOM methodpj: firstly the entangled 
discretisation and perturbative truncation errors and secondly the rather low scale at which 
continuum perturbation theory is applied. 

Phenomenology requires that scheme dependent QCD observables calculated within lat- 



tice QCD are converted to a perturbatively amenable scheme such as MS. This is often done 
using an intermediate, regularisation invariant momentum scheme (known as RI/MOM). 
The conversion of lattice results to the RI/MOM scheme is non-perturbative and the scale 
is defined by the momenta used. Obtaining precision in the lattice calculation requires 
that the momentum scale used is well below the lattice cut-off, giving the high end of the 
Rome- Southampton momentum "window" condition: 

A^cD«P^«©'. (1) 
The perturbative conversion between RI/MOM (or RI/SMOM in the case of non-exceptional 



kinematics) and MS is known for many operators to as many as three loops. Perturbative 
accuracy is dependent on satisfying the lower inequality in Eq. ([T]). The error associated 
with the low end of the window depends on logarithms the physical momentum scale, while 
the high end converges more benignly as a power of the dimensionless momentum scale {{apY 
in the case of domain wall fermions). The continuum extrapolation enabled by this work 
affords some tolerance to operating near the margins at high momentum by extrapolating 
away discretisation errors. The step scaling proposed in this work is intended to raise the 
momentum scale into a more convergent perturbative regime. 

The structure of the paper is as follows. In section [III we summarise the method of 
off-shell non-perturbative renormalisation. In section IIIII we introduce twisted boundary 
conditions for off-shell renormalisation. This is motivated in section IIII Al where we explain 
how this solves a serious theoretical problem with defining the continuum limit of operators 
renormalised using RI/MOM. The implementation of twisted boundary conditions in the 
valence sector for fermion vertex functions is discussed in section IIIIBl and we define lattice 
kinematics ( IIII CI) with which we intend to take the continuum limit ( IIII Dl) for both ex- 
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ceptional and non-exceptional momenta. Our proposed step scaling scheme, based on these 
techniques, is discussed in section HVl 

We present data for the amputated vertex functions at a single lattice spacing in sec- 
tion IV A[ In section IV CI we demonstrate practicality using existing configurations to take 
the continuum limit for the first step in the step scaling programme. 

II. BACKGROUND 

The Rome-Southampton RI/MOM approach |l| involves a simple physically defined mo- 
mentum scheme, albeit with particularly unphysical renormalisation conditions. For any 
given operator, O, a momentum configuration for the vertex function of that operator is 
chosen with momenta selected from the Fourier modes for the simulated lattice size aL^, 
where is an integer: 

27r 

, with Uf, G {0, . . .Lf,}. 

The renormalisation condition for any regularization scheme S chooses counter terms 
such that renormalised amputated vertex function, in Landau gauge and for a chosen set of 
external momenta, is equal to the tree level operator. For example, for bilinears this is 

A's''\Pi,P2) = $AUpuP2) = 0'-, (2) 

where A represents an amputated vertex function. The original paper used the amplitude 
with 

Pi=P2;g^ = 0. (3) 

The approach facilitates the conversion between different schemes because all that is required 
in any given scheme 5 is a self-consistent calculation of the relevant scattering amplitude to 
some order. By virtue of the physical definition, one scheme can be a fully non-perturbative 
lattice calculation with some lattice action at non-zero lattice spacing. Thus, within the scal- 
ing window, Eq ([T]), a lattice simulation can be matched directly to continuum perturbation 
theory without use of ill-convergent lattice perturbation theory. 

This momentum window argument motivates the neglect of discretisation errors in a 
certain limit. We shall see in section UTTI that a rigorous continuum extrapolation using the 
original method is made impossible due to the use of Fourier modes. This is because keeping 
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both the direction and magnitude of the physical momenta fixed while changing the lattice 
spacing is incompatible with the discrete momentum spectrum. 



A. Non-exceptional momentum 



In the original RI/MOM kinematic choice a single gluon can carry all the hard momentum 
through the vertex leading to infra-red effects which fall only as ^. Depending on the 
operator these effects have a non-trivial dependence on the valence quark mass that can 

nfiQ 

complicate mass extrapolations, and in some cases require pion pole subtraction l5H7|. 

RBC and UKQCD have developed the non-exceptional momentum, SMOM, kinematic 
point as a preferred matching condition js]. A significant gain comes from the extra sup- 
pression of non-perturbative effects. The perturbative expansion must be calculated for 
non-exceptional momenta; this has been performed to one loop for bilinear operators 
and for Bk |io|. Two recent publications have extended this to two loops for two different 



schemes for the quark mass 
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12| . The non-exceptional kinematic point for bilinears has 



pI=pI = {Pi-P2f- 



(4) 



RBC and UKQCD have found Sj that both momentum and mass dependence are simpli- 
fied with non-exceptional kinematics. With this momentum fiow, multiple hard gluons are 
required to create a soft sub-graph and non-perturbative infra-red effects, such as sponta- 
neous chiral symmetry breaking, fall as a higher power of the external momentum, ^ 13 |. 
As our aim is to carefully study vertex functions while reducing the volume, we require a 
momentum scheme that is not sensitive to condensate physics. 



B. Volume source NPR 



The original method used a single point source to calculate momentum space Green's 
functions. The volume source technique was developed by QCDSF [2I. The attraction here 
is to evaluate the amputated vertex function with the operator insertion averaged over all 
lattice sites. We shall see it is easily possible to obtain 0.1% statistical errors with volume 
source techniques, and even smaller if hundreds of configurations or larger volumes were 
used. With this statistical precision systematic effects like 0(4) breaking lattice artefacts 
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are visible and are in fact the dominant systematic error. These must either be included in 
the error analysis or better yet addressed using the techniques of this paper to escape the 
Fourier constraints. 

We use ^, j to represent colour indices, and a, (3 to represent spin indices. We define the 
four momentum source, used on a Landau gauge fixed configuration, as 

7]p{x) = e'P^'^'SijSafs, (5) 
where the (dimensionless) momenta are ap^ = Ufj,^. On a given gauge field U^{x) we solve 

^Ddwf(2/,a;)G(a;,p) = T]p{y), (6) 

X 

and -Ddwf is the Domain Wall fermion matrix 3|. One propagator inversion for each leg 
momentum is necessary, but this cost is more than offset by the gain in statistical accuracy 
from the volume average. 

To compute the external legs required for amputating the vertex functions we require the 
momentum space propagator 

G(pi,pi) = 5^e-^^^-G(x,pi). (7) 

X 

We form phased propagators for each momentum: 

y 

We select two independent momenta pi, and p2i and form unamputated bilinear vertex 
functions Vy for Dirac structure F: 



(9) 



^75(G"(x,pi))t75FG'(x,p2 

X 

and also for four quark operators: 

^rr'(Pi,P2) = (75(G"(x,pO)^75rG"(x,p2)),^.„^ {^,{G\x,p,)y^,TG\x,p,))^^^^^^. (10) 

X 

Here, external colour and spin indices are left free for off line amputation and projection. 
This allows us to define the amputated Green's function as, for example, 

nr(p) = (G-i(pi,Pi)Vt(pi,P2)75[G-'(P2,P2)]h5) (11) 



where p'^ = Pi = P2 = {pi — P2Y, and we take q = pi — P2- The bare vertex amphtudes are 
then obtained by projecting the amputated Green's functions onto their tree level values, 
for example 



III. CONTROLLED CONTINUUM EXTRAPOLATION 

We will motivate and introduce the use of twisted boundary conditions in the context of 
off-shell renormalisation. This will enable us to develop a framework for better controlled 
continuum extrapolation than is possible using Fourier modes. 

A. Continuum extrapolation and Fourier constraints 

In continuum Euclidean space simultaneously rotating all momenta in a scattering am- 
plitude by any 0(4) matrix must give equivalent results. However, in our discrete system 
there is only H{4) symmetry and, even in the infinite volume, only momenta related by 7r/2 
rotations and reflections are equivalent. Vertex functions will receive different O(a^) errors 
depending on the direction of the momentum relative to the lattice axes. The original studies 
did not resolve 0{4) breaking lattice artefacts due to statistical imprecision. Different 0{4) 
equivalent but ^(4) distinct momentum configurations, given by different Fourier modes, 
were treated interchangeably and simply averaged together. This is backed to some extent 
by the Rome-Southampton window argument for a region of safe operation with respect to 
lattice artefacts. 

When the method outlined above is combined with a continuum extrapolation, the naive 
use of 'H(4) inequivalent momenta as interchangeable is theoretically unsound because 0(4) 
breaking is a component of the lattice artefacts that should be removed by the extrapolation. 
If the orientation of momenta relative to lattice axes does not remain fixed as a is changed, 
an extrapolation in a is not formally valid as one uses a different observable for each lattice 
spacing (differing due to 'H(4) inequivalent momentum choices). The Fourier constraints 
make it difficult to select the same physical momentum simultaneously on more than one 
lattice spacing, and these differing 0(4) breaking lattice artefacts necessarily enter each 
data point in a continuum extrapolation of data renormalised using the Rome-Southampton 




(12) 
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method. This introduces an intrinsic systematic uncertainty at a level set by the size of 
0(4) breaking artefacts. 

0(4) breaking effects were argued to be ignorable as they were not statistically resolved. 
This justifies smooth interpolation in to obtain matched physical momenta. 

It also somewhat justifies an (otherwise dubious) continuum extrapolation. We measure 
the vertex function observable on different lattice spacings at the same physical momentum 
magnitude, but must pick 7/(4) inequivalent momenta to satisfy Fourier constraints. This 
results in incorrectly parametrised lattice artefacts, however 0(4) breaking effects can be 
safely ignored when they are substantially smaller than statistical errors. 

With the volume source technique, however, 0(4) breaking effects are well resolved and 
are a dominant systematic error. This compromises our ability to perform a continuum 
extrapolation, and these considerations represent a serious problem that we address in this 
paper with the use of twisted boundary conditions. 



B. Twisted boundary conditions 



From the discussion of the previous subsection IlII Al we see that it would be rather better 
to crisply remove lattice artefacts by continuum extrapolation, and only apply perturbative 
matching in the continuum limit. This is difficult when constrained to use Fourier modes. 
Applying twisted boundary conditions to propagator inversions enables arbitrary momenta 
to be used. The twisted boundary approach differs by only finite volume effects from a 
simulation on a much larger lattice with an exact Fourier mode of the same momentum. 

In order to simultaneously satisfy the constraint (|4]) and only use one momentum direction 
we propose choosing one kinematic satisfying (jl]) and use twisted boundary conditions 
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30|, to vary the magnitude of the momentum. The twisting technique has been used 



to insert arbitrary three-momenta in form factor calculations with the success demonstrated. 



for example, by the direct comparison of refs 
the same configurations and action. 



31 



32| to refs |30|, l33|, ISJ] • which make use of 



25|, and 



In this paper the technique is used to allow arbitrary four-momenta to be used 
hence allow to rigorously disentangle discretisation and perturbative truncation effects for 
the first time. Of course, the lattice artefacts are still present at finite lattice spacing but, up 



to finite volume effects, these become the same lattice artefact at each lattice spacing. This 
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is not the case if different Fourier modes are used. The goal is to enable us to determine the 
non-perturbative anomalous running in the RI scheme in the continuum limit prior to the 
use of continuum perturbation theory. 

We now formulate the momentum space propagator calculation with twisted boundary 
conditions for use in off-shell renormalisation. Let the quark fields satisfy twisted boundary 
conditions q{x + L) = e'^g(x) and define 28| 

q(x) = e-'^^'qix) (13) 

with aBf^i = so that q{x) satisfies periodic boundary conditions. This transformation 
changes the continuum Dirac operator: 

D = {0+m) ^ D = {0+0+m) (14) 

D has inverse G and D has inverse G; they are related by 

G{x,y) = e'^^''-y^G{x,y) (15) 

using translational invariance which we should recover after gauge averaging. Let 

G{z,p) = J2G{z,x)e^^^ (16) 

X 

then 

J2Diy,z)Giz,p) = e'Py (17) 

z 

So inverting the twisted Dirac operator ( lTi|) with a momentum source gives G{z,p). Note 
that G{z,x) satisfies periodic boundary conditions so that in the source term, e*^^, ap^ = 
is a Fourier mode. G{z,p) is related to G{z,p) via, 

G{z,p) = ^e-'^("-")G'(2,x)e'P" = e-'^'G{z,p + B). (18) 

X 

The propagator can be obtained from, 

G{p, p) = J2 e"''"G'(^, P) = Y1 e~'^''^^'^^"~"^Giz, x) = G{p + B,p + B) (19) 



z.x 



Thus the net effect of twisted boundary conditions is to shift the momentum p io p + B 
where B is arbitrary. In order to compute the non-exceptional Green's functions observe, 

^^(^1,^2) = 5^e-''P^(='-^)G(x,z)re-'^'^(^-^')G(z,2/) = Y,^,e''"'G{z,p,)^^^Te-'P^'G{z,p2) 

(20) 
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x,y,z 



Using the inverse of (fT8|) Gr{pi,P2) the vertex function is seen to be, 
x,y 

(21) 

We consider the values for pi, p2, Bi and B2 in the following subsection. 



C. Kinematics 

P 

z 

i 




Figure 1: Non-exceptional momentum configuration used in this paper. The momenta pi, p2 and 
(pi ~ P2) must be equal in magnitude and are represented as an equilateral triangle with vertices 
touching the px, Py, and pz axes. With twisted boundary conditions the sides of this triangle can be 
continuously scaled allowing both smooth interpolation of momentum dependence and controlled 
continuum extrapolation. This is very much in contrast to the situation that arises when using 
only Fourier modes. 

We select Euclidean momenta in the direction of pi = (—1, 0, 1, 0) and p2 = (0, 1,1,0) up 
to H{A) symmetry operations. This choice minimises S*^ = XliPf which we take as measure 
of discretization errors, subject to the constraint that Pi/p2/{pi — P2) be H{A) equivalent 
momenta. It is certainly possible to rotate some combination oi Px,Py,Pz into the temporal 
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direction breaking this i?(4) equivalence, and this would form an interesting possibility to 
demonstrate universality of the continuum limit in a later work. 

The non-exceptional Pi/p2/{Pi ~ P2) kinematics can be represented by an equilateral 
triangle in momentum space with vertices lying on the Px,Py,Pz coordinate axes (figured]). 
Continuous dilation of this triangle can be performed by using twisted boundary conditions 
with a fixed orientation of the pi , p2, pi — P2 triangle. This continuity allows one to pick 
a fixed MOM observable as the lattice spacing is varied. In particular, we choose parallel 
twisting vectors in the directions Bi = {—9, 0, 9, 0) corresponding to pi and B2 = (0, 9, 9, 0) 
corresponding to p2- We continuously vary 9 to move the vertices of the triangle along the 
axes. 

The non-exceptional momenta are chosen to have small discretisation errors by spreading 
the power across multiple coordinate axes. In general for any non-exceptional configuration 
there will be O(a^) errors. For exceptional momentum amplitudes we take vertex functions 
with pi on both legs. As the direction is kept fixed the vertex function is a smooth function 
of (ap)^ and interpolation in momentum is possible. Extrapolation to the continuum limit 
at fixed physical momentum is now theoretically justified. 



D. Continuum extrapolation 

Between different lattice spacings we compare only renormalised quantities and their 
ratios. This is quite natural as we desire to take the continuum limit of the product of a 
bare operator matrix element and its RI scheme renormalisation constant at some accessible 
scale p. The scaling factors for each operator to match between different /3's introduced in 



23| are then not required. The aim of the latter part of this paper will be to raise, in the 
continuum limit, this scale p to one {s^p where s is a scale factor) where there is better 
perturbative convergence. 

We eliminate Zq via the vertex function of the (conserved) vector or axial currents: 

D r N AA{p,a,m) Zo{p,a,m) 

Roip, a, m) = — = 22) 

Ao(p,a,m) Za 

and this ratio will be extrapolated to the chiral limit to produce a renormalization constant 

Zci{a,p) = Za lim Ro{p,a,m) (23) 

m— 5-0 
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where, for convenience, we use the vertex function of the local rather than conserved axial 
current and eliminate this with a previously computed Za- We then continuum extrapolate 
the product of Zq with (O): 

= lim(C)Zo. (24) 

a— >0 

This gives ((9) in the RI/MOM scheme at a scale determined by the physical momentum 
chosen which must be the same scale on all the lattices used to take the continuum limit. 

We also consider the factor required to convert (C) at scale p to scale sj) where s > 1 is 
a scale factor. This is 

So(p, sp,a) 

and its continuum limit is 

(Joij), sp) encodes the scale dependence of the RI scheme in the continuum limit expanded 
around the scale p. We can then assess the degree to which the running around our match 
scale is perturbative without risk of confusion by possible lattice artefacts. By picking a 
well defined momentum orientation for the amplitude as a function of a, we remove the 
ambiguity in taking the continuum limit that arises when selecting different Fourier modes 
for each value of /3. 

We have demonstrated how to obtain a controlled expansion in powers of the lattice spac- 
ing for amputated vertex functions and gain better control of lattice artefacts via continuum 
extrapolation. 



Ro{sp,a,m) 

hm — , 

m^o Rq[p, a,m) 



hm T^oip, sp, a). 

a-i>0 



(25) 



(26) 



IV. STEP SCALING 



QCD perturbation theory at lattice scales is not rapidly convergent, and a means to 
increase this scale without applying brute force to raise the lattice cut-off is required. Step- 

n n 

scaling [1^, 115| is the natural approach to do this where, in a series of simulations, the 
physical volume is reduced as the lattice spacing is reduced to enable the renormalisation 
scale to be raised without the cost associated with simulating all scales on a single lattice. 
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A. Step scaling background 



Step scaling has been well developed in the Schrodinger functional approach 
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201 ] ■ A finite volume scheme based on Creutz ratios has also been recently developed 
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In this paper we seek to develop a related approach based on the Landau gauge fixed 
RI/MOM method. Our approach does not tie the renormalisation scale to the volume, 
and therefore also avoids the need to fine tune bare couplings to precisely match a volume 
sensitive renormalised coupling between simulations. 

After reviewing existing step scaling methods we will attempt to address perturbative 
errors associated with the low-end of the Rome-Southampton scaling window by raising the 
scale at which we match to perturbation theory by introducing a new step scaling scheme. 
In order to use the same notational conventions as the literature on step scaling in the 
Schrodinger functional scheme, for the rest of this section let L represent the length of 
the lattice in physical units rather than units of lattice spacing. This departs from our 
convention elsewhere in the paper, which takes L as an integer lattice size and is consistent 
with notation typically used for RI/MOM renormalization. 



1. Schrodinger functional 

In step scaling there is typically a renormalisation condition imposed that is physically 
defined in a fictional finite volume universe. In the case of the Schrodinger functional this is 
the coupling g{^) at a scale = 7; defined by the volume. A second simulation, at the same 
bare coupling and with larger volume sL, then gives the coupling at a scale = f = ^ 

^(/i',a) = S(s,^(-^),-^). 

Here s > 1 is a scale factor, and is typically s = 2. The continuum limit 

(y{s,g{\)) = limS(s,^(y),y) = \img{fi' = ■^,a), 
L a~^0 L L a-^o sL 

is taken while holding the measured coupling g{fi = fixed. Hence g{^ = j^) has no 

a dependence and the trajectory to the continuum limit leaves the renormalised coupling 

g{fi = j) lattice artefact free at (and only at) the scale fi = j^. 

However, g{fi', a) is lattice spacing dependent, with discretisation errors that are removed 

by taking this continuum limit. Thus g{fi = j^) is the observable quantity used to define 
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the lattice spacing and hence physical volume. The Schodinger functional is particularly 
economical since this is in fact the most infra-red scale accessible within the simulated 
volume. 

When two consecutive scale evolution steps are considered, two sequences of simulations 
must be performed to determine the evolution of the coupling from scale fi'^ fii, and from 
f^2 ~^ in such a way that fi2 = f^i- Since different quantities g{fii) and g{fi2) are held 
fixed for the two sequences of simulations (and effectively determine the lattice spacings), 
satisfaction of the constraint fi2 = /^'i in the continuum limit is ensured by defining the 
scaling trajectory for the second sequence such that g{fi2, o) = a — 0). 

Since each simulation must correctly describe the length scales associated with the quan- 
tity used to determine trajectory to the continuum limit, carefully changing this between 
steps was a key component allowing the volume to be reduced. This is a feature we intend 
to reproduce in our method. 

In the Schrodinger functional g{j) defines the scaling trajectory and is directly coupled 
to the volume. Taking the continuum limit holding g{j^) fixed then requires fine tuning of 
(3 to exactly match the desired non-perturbative coupling - and hence to match physical 
volumes defined by using g{j^) to set the scale. This is a fine tuning step that our proposal 
below avoids. 



2. Previous work on step scaling RI/MOM 



The possibility of step scaling with RI/MOM has been previously studied 22|, 23|]. This 
work used a series of quenched configurations where the ratio of lattice spacings had been 
tuned to be precisely a factor of two to obtain aligned Fourier modes. This fine tuning is 
expensive in a dynamical simulation. Use of power counting and the Rome- Southampton 
scaling window was made, rather than the controlled continuum extrapolation at fixed phys- 
ical momentum introduced in this paper. Free parameters for each operator were introduced 
to match the renormalisation constants for different /3's corresponding to the (possibly non- 
perturbative) anomalous scaling of the operator with the lattice cut-off. This excellent start 
was not easily developed into a practical technique. 
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B. RI/MOM Step Scaling Proposal 



In order to develop a practical step scaling scheme we require several pieces. Firstly, 
we have defined continuum limit non-perturbative evolution ratios that can be used to 
access higher momentum scales. Step scaling will combine a sequence of these objects, each 
determined inexpensively compared with a brute force method. 

We believe that three technical advantages we have over the earlier attempt at RI/MOM 
step scaling make the task more tractable. The first, volume sources, reduces statistical 
errors greatly giving precision to the approach. The second, non-exceptional momenta, 
renders mass dependence and infra-red behaviour in benign; this greatly assists with 
both chiral extrapolation and matching between ensembles, in addition to finite volume 
sensitivity. These were described above. 

In this paper we have introduced a third advantage: using twisted boundary conditions 
to select arbitrary four momenta. The direction of the scattering momenta can be kept 
fixed relative to lattice axes and arbitrary values of can be chosen. This is a Good Thing 
because it allows the same physical momenta according to the lattice symmetries i?(4) to 
be chosen on each ensemble. The vertex function observable will only then have a valid 
expansion in powers of the lattice spacing a. This approach enables precise matching of 
momenta between ensembles, and precise continuum extrapolation. 

We will keep the off-shell momentum scales hard enough that the physical volume should 
not be resolved. In this way the the calculation will be finite volume effect safe, and the 
perturbative matching will use the standard infinite volume perturbation theory already 
available to high order for off-shell renormalisation. These are significant advantages over 
the Schrodinger functional. 

For the scheme to be affordable it is therefore necessary that it be able to operate in a 
small physical volume. For offshell amplitudes this is the case whenever the virtuality is too 
hard to resolve the finite volume 




(27) 




(28) 



This latter condition is likely possible to meet using modest lattice volumes ^ < 16 at all 
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stages of the calculation. We will compute vertex functions on lattices at different values of 
P with overlapping scaling windows. This will then enable the determination of a continuum 
step scaling function for some operator O at a physical scale chosen in the window Eq f l28p 
to that at a scale larger by a scale factor s. This is precisely (To{p,sp) from the previous 
section. 

These can be combined in a series of non-perturbative steps to a high scale, with pertur- 



bative conversion to MS is then applied only where it is well convergent. For example, 



(29) 



X 



1 + c^^^^^^^asiix = s» 

It appears clear that Eq (125]) can be satisfied at reasonable expense. For example 16'^ domain 
wall fermion simulations will be inexpensive with the next generation of supercomputers; 
multiple such ensembles dedicated solely to the renormalisation of lattice operators are 
quite affordable. This is particularly helped by the relatively benign mass dependence of 
non-exceptional momentum vertex functions. It is also clear that on these proposed small 
ensembles all hadronic quantities must be avoided to ensure finite volume safety of the 
analysis. For the determination of mres, the PCAC is an operator relation and so it is 
immaterial whether a finite volume or a physical pionic state is used to determine the ratio. 
Determining the lattice spacing (or more specifically the ratio of the lattice spacing to that 
of an earlier larger volume simulation) in a finite volume safe manner is discussed below. 



1. Trajectory to continuum limit 

We now define how the continuum limit of I1q{s"'~^p, s"'p, a) should be taken to obtain 
aQ{s'^~^p, s'^p). In order to avoid the fine tuning problem we seek a quantity q associated 
with a continuously variable length scale Lg significantly shorter than the lattice extent L to 
determine the trajectory to the continuum limit. We plan to reduce the volume successively 
with each step, and a different definition of the lattice spacing (e.g. a different Lg) must be 
used for each step. 
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In the Schrodinger functional method this is done by fine tuning (i.e. constraining) 
g{fin+i) for each non-zero lattice spacing entering the continuum extrapolation of cr"+^ to 
match the continuum limit of the high scale coupling obtained from the previous step, 
lima^o g{fJ''n, 0,)- We seek a the family of scale determining quantities {qn} that each be used 
to determine the lattice spacing in the corresponding n-th step of the step scaling scheme. 
These will play part of the role that g{fi = j^) plays for the Schrodinger functional. 

The continuum limit of successive quantities g„ in the sequence will be determined non- 
perturbatively as we proceed. Previously determined values will be used to define the tra- 
jectory to the continuum limit for successive steps. As a concrete example our initial results 
we will take g„ to be the static inter-quark force and so the determination of the lattice 
spacings corresponds to a family of Sommer scales based on the static potential. In the 
Schrodinger functional the quantity g„ corresponds to the coupling g{-^). 

At step n we will determine the continuum limit of Qn+i while determining the lattice 
spacing a*-^"-* from and thus must maintain the constraint g„(a) = while taking the 
continuum limit: 



so that, for example, we determine the continuum limit of the inter-quark force at a reduced 
length scale in one step, and then reuse this to constrain the trajectory to the continuum 
limit in the next step. Using this trajectory we define the scale evolution functions a". 



where is determined by the previous step. 

When g is a function of a continuous scale Lg <^ L, this distance can be varied post 
simulation and without finite volume effects, and this allows us to avoid fine tuning /3. 
Ultimately this process will become difficult as all dependence on lattice spacing becomes 
logarithmic and precision quantities are required; however, there may also be less appropriate 
but expedient choices that enable immediate progress for a few steps away from our relatively 
coarse simulations. We consider two possibilities. 




(30) 




(31) 
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Static potential 



We note the static potential has been measured successfully over a large range of length 



35|. This involved the 
. This 



scales in the quenched approximation with the Wilson gauge action 
use of a shorter length scale ~ 0.26fm than the more common Tq ~ 0.48fm 
calculation is also interesting because it successfully fits log ^ as a polynomial in /3. Such 
an approach may ultimately assist continuing lattice spacing determinations to increasingly 
fine and small volume simulations. 

We consider a sequence of scales, of the same class as the Sommer scale 

r2F(r„) = (32) 

The Sommer scale tq takes 

Co = 1.65 (33) 

Thus a step scaling scheme with scale factor s can then be defined choosing pn = s"'p and 
Tn = ^ SiS follows: 

• Determine (j{po,Pi) in continuum limit holding roPo fixed such that ro-F(ro) = Cq 

2 

• Determine Ci = ^F{^) in continuum limit holding ro fixed 

• Decrease L by ~ -1 without fine tuning 

• Determine cr(pi,p2) in continuum limit holding ripi fixed such that r^F(ri) = Ci 

2 

• Determine C2 = ^F{^) in continuum limit holding ri fixed 

• Decrease L by ~ without fine tuning 
etc... 

Following a rule of thumb that r < ^ for the static potential should ensure finite volume 
safety and enable simulation down to L ~ 0.75fm when using scales similar to to set the 
scale. Eventually this will become imprecise when entering a region where the potential 
runs logarithmically. However we believe a substantial benefit is already achievable prior to 
addressing this issue. 
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Alternative schemes 



In principle, we can use the momentum dependence of the off-shell vertex functions 
themselves to match the lattic scales between different ensembles. This may assisted by the 
improved techniques of this paper, particularly the tunable momentum scale which allows 
us to select the length scale that sets lattice spacing. This will be the subject of further 
study. 

If we accept fine tuning, it is also possible to combine step scaling of off-shell vertex 
functions with existing finite volume schemes such as the Schrodinger functional or the 
Wilson loop scheme used to match lattice spacings. 



V. RESULTS 



We use the domain wall fermion action |3i, and 16 x 32 ensembles with L., = 16 and with 
the Iwasaki gauge action {4]. We use two ensembles sets in this paper, with /3 = 2.13 and 
/3 = 2.23. The /3 = 2.13 ensembles were studied extensively in Ref. 37||. The three /3 = 2.13 
16^ ensembles used in this section have strange quark mass rrih = 0.04 and degenerate light 
quark masses m; = 0.01, 0.02, 0.03. The second ensemble set with /3 = 2.23 is not previously 



published and was used as part of a parameter search, made in the style of [38|, |39[ , prior to 
our 32'^ simulations 40| which lie nearby in parameter space. There are two ensembles with 
rrih = 0.04 and mi = 0.01,0.02. There are around 2000 trajectories in these ensembles. 



For /3 = 2.13 the extrapolations to rriq = —rrires were made with m 



res 



0.00305, Za 



0.7161(1), and the lattice spacing was taken as = 1.729 GeV 4]|. We find in this 
paper that (3 = 2.23 corresponds to an inverse lattice spacing of around ~ 2.14GeV, we 
estimate mres — 10"'^ and Za = 0.740 by interpolating between nearby values for (3. This is 
not ideal, but adequate for the purposes of this demonstration. 

All of the vertex function data was computed using the volume source method with 
twisted boundary conditions as described in section IIII Bl using 20 gauge configurations for 
each mass, see tables [T] and [TTl with momenta in table IIIII 
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arrig Range A 
0.03 800 - 2320 80 
0.02 1000 - 2520 80 
0.01 1000 - 2520 80 

Table I: /3 = 2.13 lattice giving the range of molecular dynamics time and the separation A between 
gauge configurations used in this work. Because so few configurations are needed with the volume 
averaging a large A was chosen to minimise auto-correlation effects. 

arriq Range A 
0.02 1000 - 1760 40 
0.01 1240 - 2000 40 

Table II: /? = 2.23 lattice giving the range of molecular dynamics time and the separation A 
between gauge configurations used in this work. 

A. Vertex functions with tvi^isted boundary conditions at a single lattice spacing 

In this section we focus on the data obtained using the /3 = 2.13 ensembles and display 
the quality of data obtained on a single lattice spacing. 

Figure [2] contains a comparison of the traditional volume averaged Fourier mode NPR 
to our new twisted boundary condition technique. We show the vector and axial vertex 
functions at fixed quark mass, rriq = 0.03, and with exceptional momentum kinematics on 
the /3 = 2.13 lattice, using 10 configurations. 0(4) breaking lattice artefacts produce a large 



(3 api 




ap2 




Ran^ 


;e 


2.13 {0,x,x, 


0) (- 


oc ^ ^ X ^ 


0) 


X = x(l-875) ^ 


X = f (2.75) 


2.23 (0,x,x, 


0) (- 


oc ^ ^ X ^ 


0) 


X = ^(1.5625) ^ 


X = f (3.125) 



Table III: The momenta at which we calculated the vertex functions. We note that the momentum 
^(0,1.875,1.875,0), for example, can be reached using 'base' momentum ^(0, 1, 1,0) and twist 
^(0,0.875,0.875,0) or base ^(0,2,2,0) and twist ^(0,-0.125,-0.125,0). We chose the base 
momentum such that the magnitude of the twisting component was less than ^(0.75), but such a 
choice is entirely arbitrary. 
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Figure 2: The axial (red) and vector (green) vertex functions on the /3 = 2.13 lattice at niq = 0.03, 
see Section IV Al for more details. Comparing untwisted (squares) with twisted (connected) data; 
twisting completely eliminates the 0(4) breaking scatter. 



scatter in the data using the traditional approach and this is scatter completely removed by 
the twisted boundary conditions technique. Of course, lattice artefacts are still present and 
of the same size, but the key point is that we can now vary /3 while looking at the same 
off-shell momentum in order to extrapolate these away in the continuum limit. 

For the exceptional kinematic configuration we use the projectors of js] and the pertur- 



bative running and matching calculated in 



£, 45| 



42| for mass and 43| for the tensor current. For 



We denote these exceptional momentum schemes 



Ovv+AA we use the results of 
RIMOM. 

For non-exceptional kinematics we compare the two schemes introduced in {q]. The first is 
the scheme of jsl which corresponds to choosing projectors 0^ for the vector vertex function 
and ^9^75 for the axial vertex function. We denote this scheme SMOM-^. The second uses 
7^ and 7/^75 for vector and axial vertex functions, we refer to this scheme as SMOM-7^. 
The one-loop matching and two-loop anomalous dimensions for tensor current and mass are 
given in {q]. These results have recently been extended to two loop matching and three 
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loop anomalous dimensions Q. For Ovv+aa we use several unpublished new schemes 



10|, and we thank them for their private 



and perturbative results by Sachrajda and Sturm 
communications. 

Figure |3] displays the vertex functions of all the operators analysed in this paper (bilinear 
Vector, Axial Vector, Scalar, Pseudoscalar, Tensor and four quark operator Ovv+aa) for 
the P = 2.13 ensemble, for both non-exceptional and exceptional momenta. We will return 
study each of these structures in turn and in more detail below. 

We wish to eliminate the wavefunction renormalisation by taking ratios ^ and where 
Aa is the vertex function of the local axial current, and its renormalisation is eliminated 
using Za previously determined from the ratio of matrix elements the local axial current 
and the conserved axial current From the Ward identities in the limit of small mass and 
large momentum we should find that Zy — Z^ and Zg = Zp. These are well supported 
by our data for non-exceptional momentum but, over our momentum range, not held for 
exceptional momentum. In order to compare with previous results we adopt, for exceptional 
momentum, the the strategy of eliminating the quark field renormalisation using 

|^ = ^(AA + Ay), (34) 

where factors of Za are then multiplied out using a previously calculated value. For ex- 
ceptional data this matches central value conventions with 8|. The 2% difference between 
|(A^ + Ay) and A^ is a chiral symmetry breaking uncertainty that propagates globally in 
operator renormalisation analysis as a systematic uncertainty (and is doubled for four quark 
operators). This uncertainty is eliminated for the non-exceptional momentum kinematic HI 
For the non-exceptional case we use simply = A^. In this case, we have 

ZttiZa = -T^, (35) 

— = — (36) 

Za At 

„ "^Ovv+aa ^^A (o>7\ 
Zbk = ^2 - 7 • 137 j 

^A ^^Ovv+AA 

Our results for RIMOM, SMOM-^, and SMOM-7^ schemes for the wave function renor- 
malisation determined from the axial current vertex function (non-exceptional) and average 
of axial and vector vertex functions (exceptional) are displayed in figure O Here, the running 
at this lattice spacing is poorly described by continuum perturbation theory, and is asso- 
ciated with the momentum dependence of the exponent for binding of light modes in the 
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Figure 3: All the vertex functions that have been calculated on the /3 = 2.13 lattice. Non- 
exceptional and exceptional vertex functions in the chiral limit. Note the great reduction in chiral 
symmetry breaking effects (A^ — Ay and Ag — Ap) at non-exceptional momentum. 



22 



1.16 



1.14- 




1.12 



(ap) 



(a)A^, Ay exceptional 



1.06 





1.045 




(b)Ayi, Ay non-exceptional 7^ scheme 



Figure 4: Zoomed in view of Ka and Ay in the chiral limit with non-exceptional and exceptional 
kinematics from figure [3l Note the expanded scale in the non-exceptional plot. The roughly 2% 
difference in the exceptional case is eliminated, even at the lowest momenta, for non-exceptional 
kinematics. 
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Figure 5: Zq on the /? = 2.13 lattice in a MOM scheme (red) and with the perturbative running 
divided out and converted to MS (blue). Zg contains strong lattice artefacts for DWF and is 
opposite to the continuum running. We find later that this becomes consistent in the continuum 
limit. 
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46| . We note that this Zq is cancelled in the 



fifth dimension in the domain wall formulation 
ratios above when treating other operators. Naturally, one expects that such discretisation 
effects will be removed if a controlled continuum limit is taken, and this will be revisited in 
later sections. 

Figures El H and |H] display the corresponding data for the mass, tensor current and four 
quark operator Oyyj^AA relevant for Bk- As promised, the data is extremely precise and 
0(4) breaking discretisation effects do not introduce scatter in the data as the momentum 
is changed because we are always selecting the same lattice momentum with our twisted 
boundary conditions. Errors are at the 0.1% level even on this relatively small volume and 
with only twenty configurations. 

In contrast to Zg, the running of Z^, and Zyy^^A is reasonably well described by 
continuum perturbation theory for larger values of p^, even on this non-zero lattice spacing. 



This is indicated by reduction of the slope of the data after perturbative conversion to MS 
(blue) especially at higher momenta. Residual j?' dependence remains after this conversion 
at the few percent scale. This refiects an admixture of lattice artefacts, and perturbative 
truncation error. The methods introduced in this paper enable the continuum limit of the 
scale evolution to be determined and disentangle these sources of error. 

Tablel^shows ou. .suits co,„pa.ed with the .suits of Ref fl aft. foUowiug the sa,ne 



procedure of converting to the MS scheme at 2 GeV and then extrapolating to = 0. 
In all cases the extrapolation was performed using a linear fit of Ao as a function of (ap)^ 
in the range 2 < (ap)^ < 3.2. This procedure is fiawed. We note that this is in danger 
of extrapolating perturbative errors into the infrared. In the next section we show instead 
how to take the continuum limit at non-zero momentum to eliminate discretisation errors. 
These numbers are thus for comparison only, and constitute a check and a demonstration 
of the precision that is possible using the techniques of this paper. The numbers do not 
represent an attempt to improve on previous estimates of renormalization constants for the 
RBC-UKQCD domain wall programme. The physical RBC-UKQCD predictions will be 
updated using the methods of this paper techniques in later works, where we will also use 
step scaling to raise the momentum scale to minimise perturbative error. 
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Figure 6: Z„i on the /3 = 2.13 lattice in a MOM scheme (red) and with the perturbative running 
divided out and converted to MS (blue). 
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Figure 7: on the (3 = 2.13 lattice in a MOM scheme (red) and with the perturbative running 
divided out and converted to MS (blue). 
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Figure 8: Zyv+AA on the /? = 2.13 lattice in a MOM scheme (red) and with the perturbative 
running divided out and converted to M S (blue) . 
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Zo Ref: RIMOM (extrap) RIMOM SMOM-^ SMOM-7^ 

Zq 0.7726(48 + 150) 0.7753(15) 0.79605(63) 0.8035(10) 0.77744(31) 
Zra 1.656(30 + 83) 1.483(28) 1.546(15) 1.5073(26) 1.5405(18) 
Zt 0.7950(34 + 150) 0.7962(13) 0.80339(62) 0.80626(73) 0.77768(20) 



Table IV: The quark field, mass and tensor current renormalization constants in MS at 2{GeV). 
Error in tlie first column is (stat + sys) error in the other columns statistical only. Note this 



work used 20 configurations at each mass whereas 



used 300 point source measurements on 75 



configurations. The third column uses the value extrapolated to zero for comparison with the 



results of 



. The others use simple interpolation to obtain the value at = 2GeV 



Ref: [8j RIMOM (extrap) RIMOM SMOM-(^,^) SM0M-(7/„7^) SM0M-(^,7^) SM0M-(7^,^) 
0.9276(52 + 220) 0.93330(73) 0.92994(54) 0.97737(59) 0.93406(60) 1.0233(21) 0.8903(19) 

Table V: BK renormalization constant with same parameters as above. The extrapolated RIMOM 
is for comparison with 8] the rest of the measurements use the interpolated value at p'^ = 2GeV. 

B. Scale determination in a small volume 

Determining the lattice spacing for our /3 = 2.23 ensemble is a useful test case for the 
methods of section IIVB II As a preliminary we have calculated rc for our 16^ lattices. We 
compute timelike Wilson loops with four hits of APE smearing, smearing parameter 2.5, in 



35| is here required for t 



le Iwasaki 



the spatial direction. The tree level improvement of 
gluon action. We compute the tree level Wilson loop using code developed in [47[ to obtain 
Table NH where the potential includes the self energy part. 

aV{Ra) = - lim 1 "^^^^^^ = CFV{R)g^ + 0{g^) (38) 

By demanding 

F{ri) = {V{r) - V{r - a))/a = F^ree = ^ (39) 

we solve for rj using Table IVll This approach leads to much reduced lattice O(a^) errors 



35l | |48|. Further, using the force instead of the potential removes the self energy part and 
reduces the linear term, err, to a constant. Fitting the Cornell potential to find rc requires 
a. However, a is dependant on the long distance behaviour of V{r), and to constrain it one 
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R 



Vw 



Vi 



(1,0,0) 


0.166667 0.08963(1) 


(2,0,0) 


0.209842 0.12569(2) 


(3,0,0) 


0.225186 0.13879(3) 


(4,0,0) 


0.232442 0.14506(4) 


(5,0,0) 


0.236630 0.14878(4) 


(6,0,0) 


0.239366 0.15133(5) 


(7,0,0) 


0.241300 0.15318(5) 


(8,0,0) 


0.242742 0.15456(5) 



Table VI: Vw is the static potential tree level part with the Wilson gluon action, Vi uses the 
Iwasaki gluon action. As R tends to infinity the difference between successive terms for Wilson 
and Iwasaki actions (the force) is the same since the self energy part cancels and the two actions 
reproduce the same IR physics. 



has to sample large r. When fitting to the force, can be obtained without including large 
distance data in the procedure, and rc computed this way is a finite volume safe observable 
with which to determine the lattice spacing. 

In order to extract V{r) from the average Wilson loops (W{r, t)) we first plot the effective 
potential 



log I 



(iy(r,t+l)). 



(40) 



{W{r,t)) ' 

as a function of t. The largest value of r used in this analysis is 4a. Excited state contamina- 
tion and statistical noise are problematic, and we use the 'black box' method 49| to define 
an improved effective mass that takes account of the first excited state. This allows us to 
extend tmin to t = 3 which gives a significant error reduction compared to the usual effective 
mass. On the /3 = 2.13 lattice we use 600 configurations for each mass, each configuration 
rotated to use 4 different time directions. For (5 = 2.23 we use 750 configurations per mass 
and again rotate the time direction. 

The equation r'^F{r) = C has solutions in the range of our data for C G (0.6, 2.0). Lower 
values give more accurate rc but these are likely to have large discretisation errors. For any 
C however should be constant, we plot this ratio as a function of C in figure fTOl and 
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determine C > 1.4 is appropriate and gives the value for the ratio of lattice scales, 



Ri 



(02, 



23) 



102.13 J 



0.652(21) 



This implies a lattice spacing of a ^ ~ 2.14 GeV for the f3 = 2.23 ensembles. 



(41) 



> 1.2 




Figure 9: Effective mass plots for (left to right) m = 0.03, m = 0.02, m = 0.01 on the (3 = 2.13 
lattice at r = 4a. Datapoints are computed via equation 1401 The red constant is obtained from a 
fit over the last three datapoints. The blue is from the black box method with tmin = 3. 
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Figure 10: The ratio '''^'"^■^^ as a function of C. The red data is computed using the effective mass 
liniimm = 4 and the blue uses the black box method with tmin = 3, and gives the same value but 
with significantly smaller error. 

We propose below a scale factor s = 1.5, and believe that at least two iterations of step 
scaling should be possible based only on the static potential. However, the static potential 
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displays percent scale errors even with many configurations and finding a more precise 
alternative would certainly be good in any case. 



C. Continuum limit scale evolution functions 

Following equation [26] we can compute cto{p, sp) using two lattice spacings. We consider 
the evolution of renormalisation constants Zg, Z^, Zx, and Zbj^ in both the exceptional and 
the non-exceptional kinematic schemes. We choose p ~ 2{GeV) and compute (To{Pi sp) as a 
function of s. 

Because of the very smooth dependence of the vertex functions computed with twisted 
boundary conditions we can perform a simple interpolation to match values of p^ on each 
lattice. One could in principle simulate at the same physical momentum on each lattice by 
choosing the twisting angles appropriately, however since the interpolation introduced only 
a very small uncertainty this was not necessary. 

For each operator we can now evaluate T^oiP = 2GeV, sp, = 1.729GeV) and = 
2GeV, sp, a"^ = 2.14GeV). Linear extrapolation in with only two datapoints is not as 
robust as one might like, but certainly suffices to demonstrate the method. Using this, we 
can obtain the continuum limit step scaling function o"o(p = 2GeV, sp) for the quark field 
(figure [TTl) . mass (figure [T^ . tensor current (figure US]), and the four quark operator Oyv+AA 
(figure [T^ . With an appropriate third and smaller lattice spacing using a correspondingly 
smaller volume we could similarly determine the next step, evolving from 3 GeV to 4.5 GeV 
and so on. 

Figure [TT] is especially important, because Zg suffers from the greatest lattice artefacts for 
Domain Wall Fermions, and is carefully eliminated in the NPR analysis of other operators. 

In the 7^ scheme, with both exceptional and non-exceptional momenta, the finite lattice 
spacing running of Zg is in the opposite direction to the perturbative prediction, however the 
running behaviour recovered in the continuum limit is close to perturbative. The determina- 
tion of Zq with domain wall fermions displays the momentum dependence of the exponent. 



a{p), of localisation in the fifth dimension 46|. As shown here, in a number of schemes, as 



long as an unambiguous continuum limit can be defined, the Domain- Wall action will then 
produce the continuum scaling behaviour of Zg and the other renormalization constants. 
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Figure 11: Quark field renormalization running from ~ 2 GeV to ~ 3.0 GeV in the three different 

schemes compared to the perturbative running in each scheme. 
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Figure 12: Quark mass renormalization running from ~ 2 GeV to ~ 3.0 GeV in the three different 
schemes compared to the perturbative running in each scheme. 
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Figure 13: Tensor current renormalization running from ~ 2 GeV to ~ 3.0 GeV in the three 
different schemes compared to the perturbative running in each scheme. 
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Figure 14: Quark field renormalization running from ~ 2 GeV to ~ 3.0 GeV in the three different 
schemes compared to the perturbative running in each scheme. 
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VI. CONCLUSIONS 



We have introduced the use of twisted boundary conditions in off-shell renormalisation. 
This enables controlled continuum extrapolation of the relevant amplitudes and rigorous 
disentangling of lattice artefacts from perturbative truncation errors in the method. 

We outlined a step scaling approach based on the scheme that will allow us to raise 
the scales at which perturbation theory is applied. The method has been demonstrated by 
taking the continuum limit of the first step scaling function in the process. We can certainly 
raise this scale from around 2 GeV to at least around 5 GeV, and perhaps higher. The 
5GeV upper scale is limited only by the range over which our current matching based on 
the static potential is likely to be precise. This is, in any case, the scale above which one 
should consider charm and bottom quark effects. Were an appropriately precise method 
for matching lattice scales (or renormalised couplings) available, the method could then be 
applied to evolve the three flavour theory to very high momentum scales, of order Mz- 

The lattice spacing matching strategy will be the subject of further study, however, 
raising the scale at which, typically two or three loop, perturbation theory is applied for 
Rome-Southampton renormalisation for much important lattice phenomenology is already a 
significant step. For example, a naive estimate of an truncation error with 0(1) coefficient 
is reduced from 3% to under 1% by raising this scale to 5 GeV. 

We plan to produce continuum limit scale evolution functions spanning the complete 
range of lattice operators covering the region 2-5 GeV. This includes all lattice bilinear and 
four quark operators, bilinear operators with covariant derivatives for structure functions 
and distribution amplitudes, and three quark operators relevant for proton decay matrix 
elements. 

Obtaining the scaling functions in the continuum limit will enable calculations with any 
lattice action to raise the reference scale at which operators are quoted in MS from 2 GeV 
to 5 GeV with substantial reduction in systematic errors. 
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